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METHOD FOR GENERATING ROCKET MOTOR DATA PROFILES FOR FIRINGS AT 

VARIOUS TEMPERATURES 

STATEMENT OF GOVERNMENT INTEREST . 
The invention described herein may be manufactured and used by or for the Government 
of the United States of America for governmental purposes without payment of any royalties 
thereon or therefor. 

BACKGROUND 

a. Field of Invention: 

The present invention relates to generating rocket motor burnback data, pressure versus 
time, and thrust versus time data for various firing temperatures, using generally available 
specifications of the rocket motor and test firing data. 

b. Background 

In order to simulate the ballistics performance of a particular rocket motor, it is necessary 

to have a burnback profile. A burnback profile is developed by plotting the burn surface area as 

a function of the linear web burned. FIG's 7, 9, 1 1 and 13 are examples of burnback profiles. 

The area under the curve equals the volume of propellant burned. The curve is unique for each 

different rocket motor and is a function of the geometry of the propellant grain. The burnback 

profiles are generally available only to those personnel intimately familiar with the design and 

development of the motor. Even when such detailed information is available, correction factors 

which are often unique for a given temperature must be applied to produce a valid model of 
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ballistics performance. Test data from actual firings can be used to plot the curve and produce a 
burnback profile but rocket motors are normally tested at extreme temperatures only and actual 
firing test data is generally not available for temperatures within the normal range of use, for the 
rocket motor. It would be expensive and impractical to conduct a series of test firings, for each 
rocket motor, in order to obtain representative data for a spectrum of temperatures, so that curves 
could be generated from the data, for representative firing temperatures. 

Computer models have been developed to take advantage of the available test data, for 
firings at extreme temperatures, and to produce a predicted curve. U.S. Patent No.6,290,505 to 
Garcia uses a thrust versus time curve from a test firing and proportionally expands or contracts 
the base line to an interpolated ending point for a temperature, for which the data is to be 
modeled. A computer program produces a thrust-time curve having the same basic shape and 
equivalent impulse (area under the curve). The interpolation approach used in the Garcia '505 
Patent is more of a graphical re-proportioning of the thrust-time curve. 

There remains a need for a method and a computer model which can mathematically 
develop an accurate curve and data profiles for a temperature different from those temperatures 
for which firing test data is available. 

SUMMARY OF THE INVENTION 
It is an object of the present invention to provide a method and implementation by 
computer model, which can mathematically develop an accurate curve for a temperature 
different from those temperatures for which firing test data is available. 
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It is a further object of the present invention to provide the computer model in the form of 
software which requires only readily available data from the rocket motor specifications and test 
firing data. 

In accordance with the objects, the present invention is a method for predicting rocket 
motor ballistic performance at a user-selected firing temperature. The method generally involves 
processing the pressure and thrust data through a number of iterative steps to generate burnback 
profile data. First the method uses measured pressure data and varies the coefficients of a throat 
erosion equation so that the calculated final throat diameter matches user input data. Then the 
method takes the measured pressure data and varies the reference burn rate until the final 
calculated web burned matches a user input value. Finally, the method uses the measured 
pressure data, along with the calculated burn rate and throat area data, and varies the 
characteristic velocity until the area under the burn surface versus web burned curve equals the 
user input propellant weight. The result is a burnback profile with the correct final web that 
integrates to the correct propellant weight. This burnback profile, throat area profile, and thrust 
coefficient profile are useful in and of themselves for many engineering analyses. Moreover, 
once the input burnback profile has been generated, the temperature dependent variables are 
adjusted and iterative processes similar to those used to generate the burnback profile are used to 
generate output pressure and thrust versus time data at the output temperature. A more detailed 
outline of the foregoing method steps is as follows: 

First, the method develops data profiles at a test firing temperature, comprising the 
following steps: 
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1 st : Receive input of data from an actual test firing, including firing temperature, initial 
throat diameter, final throat diameter, pressure versus time and thrust versus time and input of 
data from rocket motor specifications including propellant weight, propellant density data for a 
range of temperatures, web thickness, pressure exponent data for a range of temperatures, 
5 thermal expansion coefficient, and temperature sensitivity coefficient; 

2 nd : Process the pressure data in an equation for the change in the throat diameter, radius, 
or area, attributable to erosion by making repeated calculations, in an iterative process varying 
the coefficient in the expression for the rate of erosion, and solving for the incremental change in 
throat dimension as a function of chamber pressure during an incremental time period. The 
10 calculations are made for values of throat dimension from the known initial dimension to the 
known final dimension. The result of the process yields a selected correct coefficient in the 
expression for throat erosion as a function of chamber pressure. 

3 rd : Calculate an array of time incremental values of throat dimension with the selected 
coefficient of throat erosion. 
1 5 4 th : Process the pressure and pressure exponent data in an iterative operation, using an 

equation for the burn rate as a function of incremental chamber pressure, reference pressure, a 
reference burn rate and a pressure exponent and a second equation which gives the linear web 
burned as a function of the burn rate and the time elapsed, said process varying the reference 

burn rate, in the equation for linear web burned to solve for the total linear web burned and 

-- 

20 selecting the reference bum rate which produces a total linear web burned matching the user 
input value for web thickness. Solve the equation for linear web burned for incremental times, 
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using the selected reference burn rate to yield an array of linear web burned values at the test 
firing temperature. 

5 th : Process the pressure data with the calculated throat dimension data in an iterative 
operation, varying the characteristic velocity, in an expression for incremental propellant weight 
burned as a function of pressure, characteristic velocity and throat dimension, to solve for a total 
propellant weight burned and selecting the characteristic velocity which produces a total 
propellant weight burned matching the user input value for propellant weight. Solve an equation 
for incremental burn surface using the selected characteristic velocity and reference burn rate to 
yield an array of incremental burn surface values at the test firing temperature and solve an 
equation for thrust coefficient as a function of thrust, pressure and throat dimension to yield an 
array of incremental thrust coefficient values, at the test firing temperature. 

Second, the method develops predicted data profiles for a firing temperature selected by 
the user, comprising the additional steps as follows: 

6 th : Process the test firing temperature and the temperature selected by the user in an 
expression for linear web burned, at the temperature selected by the user, as a function of the 
linear web burned, at the test firing temperature, the difference in the temperatures and a thermal 
expansion coefficient to solve for the linear web burned, at the temperature selected by the user 
for an array of integer numbered increments. 

7 : Process the array of burn surface values, at the test firing temperature, in an 
expression for the burn surface at the temperature selected by the user as a function of the bum 
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surface, at the test firing temperature, the thermal expansion coefficient, the difference in 
temperatures and a burn surface modifier, which is given an initial default value of "1 and 
solve for an array of burn surface values, at the temperature selected by the user, for integer 
numbered increments. 

8 th : Process the reference burn rate at the test firing temperature in an expression for the 
temperature dependency of the burn rate, and solve for an array of burn rate values at the 
temperature selected by the user, for integer numbered increments. The method uses the 
pressure exponent for the temperature selected by the user and the pressure value for a next 
preceding increment, in each integer numbered increment. The initial pressure prior to the first 
increment is given a default base line value. 

9 th : Process the array of burn surface values, the array of burn rate values and the array of 
throat dimension values, in an expression for the pressure as a function of the burn surface, burn 
rate, propellant density, throat dimension, and characteristic velocity to solve for an array of 
pressure values at integer numbered increments. Each value in the array provides a pressure 
value for the next following increment of the 8 th step. 

10 th : Process the array of thrust coefficient values at the test firing temperature, in an 
expression for temperature dependence of the thrust coefficient to solve for an array of thrust 
coefficient values at the temperature selected by the user. 

1 1 th : Process the arrays of thrust coefficient, pressure and throat dimension values, in an 
expression for thrust as a function of thrust coefficient, pressure and throat dimension to solve 
for an array of thrust values for integer numbered increments. 
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12 : Process the arrays of linear web burned and burn rate at the temperature selected by 
the user in an expression for the time increment as a function of the linear web burned and burn 
rate to solve for an array of integer numbered time increments and assign the resulting time 
values to the integer numbered increments in the arrays for the temperature selected by the user. 
5 1 3 th : Repeat the steps 6 th through 12 th , in an iterative operation, varying the burn surface 

modifier, in an expression for incremental propellant weight burned as a function of incremental 
burn surface and incremental burn rate, both at the temperature selected by the user, solving for 
the total propellant weight burned and selecting the burn surface modifier which produces a total 
propellant weight burned matching the user input value for propellant weight. The method 

10 rebuilds the arrays, by repeating steps 6 th through 12 th , using the selected bum surface modifier. 
The foregoing method eliminates the need for knowing the actual geometry of the 
propellant grain and makes use of test firing data, such as, pressure and thrust data versus time, 
and measured throat area, along with general propellant grain and motor design data, such as 
propellant weight and density, web thickness, propellant pressure coefficient, thermal expansion 

1 5 coefficient, and temperature sensitivity coefficient to develop a burnback profile. Additionally, 
the method uses the user input information together with the profiles it has generated to build 
data profiles for temperatures, for which no test firing data is available. The type of input data 
required is that which is readily available to personnel involved in the analysis on a given rocket 
motor. 

20 

BRIEF DESCRIPTION OF THE DRAWINGS 
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FIG. 1 is an exemplary profile of the thrust versus time for an existing propellant rocket. 

FIG. 2 (A-E) is the computer code showing a portion of the main routine of processing to 
subroutines and presentation. 

FIG. 3(A-B) is the computer code subroutine for the iterative process to calculate the 
coefficient in the equation for the rate of throat erosion. 

FIG. 4(A-B) is the computer code subroutine for the iterative process to calculate the 
reference burn rate and for producing the web burned array. 

FIG. 5(A-B) is the computer code subroutine for the iterative process to calculate the 
characteristic velocity and for producing the burn surface and thrust coefficient arrays. 

FIG. 6(A-B) is the computer code subroutine for the iterative process to calculate the 
burn surface modifier and for producing the burnback profile. 

FIG. 7 is a burnback profile generated by the present invention for two MK 106 rocket 

motors, 

FIG 8 is a graph of thrust versus time for MK 106 rocket motors showing test firing 
results and modeled predictions using the present invention. 

FIG 9 is a burnback profile generated by the present invention for three MK 106 rocket 

motors. 

FIG 10 is a graph of thrust versus time for MK 106 rocket motors showing test firing 
results and modeled predictions using the present invention and including a test firing at the 
same temperature used by the model. 
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FIG 1 1 is a burnback profile generated by the present invention for MK 1 1 1 rocket 

motors. 

FIG 12 is a graph of thrust versus time for MK 1 1 1 rocket motors showing test firing 
results and modeled predictions using the present invention. 
5 FIG 13 is a burnback profile generated by the present invention for MK 75 rocket motors. 

FIG 14 is a graph of thrust versus time for MK 75 rocket motors showing test firing 
results and modeled predictions using the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
10 The present invention is a method for using a computer model for generating rocket 

motor profiles as in FIG. 1 , inclusive of throat area profiles and thrust coefficient profiles, all 
from readily available rocket motor specification data and test firing data. The present method 
facilitates the generation of burnback profiles and pressure and thrust versus time data at firing 
temperatures for which no test firing data is available. This data is valuable for predicting the 
1 5 ballistics performance of particular rocket motors at various temperatures. 

Typically a user who is engaged in analysis of rocket motors has access to specifications 
including web thickness, propellant weight, propellant density, initial throat diameter, thermal 
expansion coefficient, temperature sensitivity coefficient and pressure exponent. The user also 
has access to the data from one or more test firings, usually at extreme temperatures. This data 
20 comprises values for thrust at time increments and chamber pressure at time increments. 
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The preferred embodiment of the present method is implemented in software form as a 
computer program comprising a main routine and a number of subroutines, the subroutines each 
performing an algorithmic calculation for collectively generating rocket motor burnback profiles 
and pressure and thrust versus time data from the existing data. 

FIGs. 2A-2E are printouts of the main routine computer code, which begins by importing 
a data file containing the above described data: web thickness, propellant weight, propellant 
density data, initial throat diameter, thermal expansion coefficient, temperature sensitivity 
coefficient and pressure exponent data, and test firing data. 

As seen in FIG. 2A, the code verifies the presence of required data elements and calls 
upon the subroutine calculations of the present invention. 

The method then proceeds to a first subroutine which is designated as Private Sub 
SetErosion that produces a throat area profile at the test firing temperature. This first subroutine 
accounts for the fact that thrust is influenced by chamber pressure and the area of the throat, 
through which the products of combustion exit. The area of the throat changes during the 
operation of the rocket motor due to erosion. The rate of throat erosion is given as a rate of 
change in the radius of the throat as follows: 




[Equation #1] 



Where: 
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5 



Ej = Incremental radial erosion rate 
Pi = Incremental chamber pressure 
a = Empirically derived coefficient 
b = Empirically derived coefficient 

The coefficient b is set at a constant value of 0.9240 which has been determined from 
experimental tests of rocket motors. The model varies the coefficient a; however, the difference 
in the throat area profile is not significant regardless of which coefficient is varied. The throat 
radius, at a given increment of time is given by the equation: 

10 

[Equation #2] 



n = + aPi (U- 



Equation #2 incorporates the expression for the rate of change in throat radius due to 
throat erosion and yields the throat radius at a given increment of time and, when solved 
1 5 repeatedly for time increments over the period during which the motor operates, yields an array 
of incremental throat radius values. 

At a given time increment, nr 2 yields the throat area at the time increment U . The 
incremental throat area is given by the equation: 



At i = TTTi 2 



[Equation #3] 
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As shown in the relevant portions of code in FIGs. 3A and3B 5 the subroutine designated 
as Private Sub SetErosion varies the coefficient a over a range of values and calculates the array 
of incremental throat radii, repeating the process until the final throat radius matches the final 
5 throat radius input by the user. In addition to the value for throat radius, the program can be 
designed to work with other dimensions of the throat, such as diameter or area. The Private Sub 
SetErosion subroutine uses the expression for the incremental throat radius, given by equation 
#2, in equation #3, including the selected value for a and solves equation #3 for each time 
increment to produce a throat area profile at the test firing temperature. 
10 The variables set forth in the Private Sub SetErosion subroutine of FIGs 3 A and 3B are as 

follows: 

Incremental chamber pressure at time increment "f - Pressln (i) 

Time increment "f - Timeln (i) 

Incremental throat area - Atln (i) and This Value 
15 Given the foregoing, the main program code then calls a second subroutine as shown in 

FIGs. 4A and 4B and designated Sub SetBurnRate, which varies the reference bum rate in an 
iterative process until the total linear web burned at the end of the firing matches the user input 
value for web thickness. 

The incremental linear rate at which the web is burned is given by the equation: 

20 
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Pi 

Yref) 



\ n 



[Equation #4] 



Pre/ 



n 



Where: 



Burn rate at reference pressure and motor conditioning temperature 
Reference pressure (model uses 1500 pounds per square inch) 



Pressure exponent input by the user 



The total linear web burned as of time increment "?' is given by the equation: 

10 



wbi = wbi - i + 



V 



(ti~ ti- l) 



[Equation #5] 



Where: wbj = 



Total linear web burned as of time increment 



The Private Sub SetBurnRate subroutine of FIGs. 4A & 4B uses the expression for the 
incremental linear burn rate given by equation #4 and solves equation #5 repeatedly, varying the 
reference burn rate in an iterative process until the total linear web burned at the end of the firing 
matches the user input value for web thickness. The selected reference burn rate is the correct 
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reference burn rate for the rocket motor fired at the test temperature. The model uses the 
expression for the incremental burn rate, given by equation #4, including the selected reference 
burn rate and solves equation #5, for each time increment to produce a web burned array for each 
time increment. 

The variables set forth in the foregoing equations are indicated in the code of FIGs. 4A & 
4B as follows: 

Burn rate at reference pressure and motor conditioning temperature - rref 
Reference pressure - PREF 
Pressure exponent - nexp 

Total linear web burned at time increment "f - Wburnln(i) and ThisValue 
Incremental linear burn rate - brate 

Given the foregoing, the main program code then calls a third subroutine as shown in 
FIGs. 5A & 5B and designated Private Sub SetCStar which calculates the mass flow values and 
varies the characteristic velocity in an iterative process until the propellant weight burned at the 
end of the firing matches the actual value input by the user for the propellant weight. 

The mass flow at a time increment "f " is a function of the incremental chamber pressure 
and the incremental throat area. The mass flow at a time increment "f is given by the equation: 




[Equation #6] 



~ c 
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Where: C* Characteristic velocity 

g Gravitational constant 

5 

The burn surface at time increment "f ' is given by the equation: 

[Equation #7] 

Si = ir 

j8p 

10 

Where: p is propellant density at firing temperature 

The method assumes a default value of 0.065 for the propellant density although the user 
can input another value for the propellant density in the original data set. 
1 5 The propellant weight burned as of a time increment "f ' is given by the equation: 



„, „, (*-l+J*Y [Equation*^ 

Pbi = Pbi - 1 + Jo -ti-i) 
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The incremental thrust coefficient is given by the equation: 



[Equation #9] 

Qi = 

PiAti 

5 

Where: pi is thrust at increment "f 



The Private Sub SetCStar subroutine of FIGs. 5 A and 5B uses equation #6 to provide the 
incremental mass flow values, with the characteristic velocity given an initial default value of 

10 60,000. The model solves equations #4, #7, #8 and #9 for each time increment, repeatedly, 

varying the characteristic velocity in an iterative process and selecting the value of characteristic 
velocity which gives the propellant weight burned, given by Equation #8, at the end of the firing, 
matching the actual value input by the user for the propellant weight. The same subroutine also 
solves equations #7 and #9, for each time increment, using the selected value for characteristic 

15 velocity, to produce an array of burn surface and thrust coefficients, respectively for each time 
increment. 

The variables set forth in the third subroutine Private Sub SetCStar at FIGs. 5A and 5B 
are represented as follows: 

Characteristic velocity - CS 
20 Burn surface at time increment "f - Bsurfln(i) 

16 



PATENT 



83,669 



Mass flow at time increment "f ' - mpdot 

Incremental thrust coefficient - cfln (i) 

Propellant weight burned as time increment "f 5 - This Value 

The total web burned and the burn surface arrays may be combined in a graph of points at 
matching time increments to produce a burnback profile which integrates to the correct total 
propellant volume burned during the test firing. The arrays of thrust coefficient and burn rate 
data form valuable profiles. This collection of information is valuable in analysis of the rocket 
motor and also forms the basis for the prediction of the ballistics performance of the rocket 
motor at temperatures for which no test firing data is available. The method uses the 
terminology "output temperature" to refer to the temperature to be modeled as selected by the 
user and the terminology "input temperature" to refer to the temperature for which test firing 
data is available. 

Given the foregoing, the main program code then calls a fourth subroutine as shown in 
FIGs. 6A and 6B and designated Private Function GenerateOutput which calculates a prediction 
of the ballistics performance of the rocket motor at any given temperature, including those at 
which no test firing data is available. 

The web at the output temperature differs from the web, at the input temperature due to 
thermal expansion. The web, at the output temperature is a function of the web, at the input 
temperature, multiplied by a thermal expansion coefficient, a\ proportioned by the difference in - 
temperatures, as follows: 

17 
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wbout'WbtnQ + ttTout-Ttfi 



[Equation #10] 



The user has the option of selecting the web, at the output temperature, in which case, the 
thermal expansion coefficient is derived by the model using the expression: 

wb c 



out _2 



in 



Tout Tin 



The burn surface at the output temperature is given by the equation: 



Sou,, ' (' + 2 x "'from " TjSin, X BS Modifier ^ ion « « 



The incremental burn rate is temperature dependent and, for the output temperature, is 
calculated, using Equation #4, with the temperature dependent variables adjusted, as follows 



&iti = &fout 



(Pi -A 



[Equation #12] 



V Pr,// 



Where: now is the pressure exponent at the output temperature. 
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&fout = &fm x e 



Op (Tout- Tin) 



The temperature sensitivity coefficient, q> , is input by the user but has a default value 
ofO.001. 

The Private Function GenerateOutput subroutine of FIGs. 6 A and 6B defines variables as 
follows: 

Reference burn rate at the output temperature - OutBurnRate 
Reference burn rate at the input temperature - InBurnRate 



Temperature sensitivity coefficient - SigmaP 
Output temperature - OutTemp 
Input temperature - InTemp 

The values for characteristic velocity and the rate of throat erosion are applied to the 
equations which determine the relationship of the variables, using the values for the output 
temperature to predict the ballistics performance of the rocket motor at the output temperature. 
Equation #7 can be expressed as follows: 



e-E 




[Equation #13] 



Where: 



Propellant density at the output temperature 
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Normally, the final-throat area is not temperature dependent; however, the method allows 
the user to select a different final throat area, in which case the method provides an expression 
for the incremental throat area to be substituted for the incremental throat area taken from the 
throat area profile generated from equation #3, as follows: 

5 

(Atim- Ato) 

Atouti= Ato\ — —(AtF- Ato) 

(AtF- Ato) J r 

v 7 [Equation #14] 

Where: 

Ato = Initial throat diameter 

Atim = Incremental input throat diameter calculated from input data 

AtF = Actual final throat diameter 

AtFmt = Output final throat diameter 

The expression for the mass flow, given by equation #13, is substituted into equation #6, 
which is rewritten as follows: 



Pi = 



v Atouug J 



[Equation #15] 
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The variables set forth in the Private Function GenerateOutput subroutine of FIGs. 6 A 
and 6B, are as follows: 

Linear web burned at increment "f , at output temperature - WburnOut (i) 
Web depth at input temperature - InWebDepth 
Web depth at output temperature - OutWebDepth 
Incremental burn surface at output temperature - BsurfOut (i) 
Incremental chamber pressure at output temperature - PressOut (i) 
Incremental throat area at output temperature - AtOut (i) 
Initial throat area - AtO 
Final throat area - AtF 

Final throat area at output temperature - AtFOut 

The thrust, at the output temperature, depends on the thrust coefficient according to 
Equation #9. The temperature dependence of the thrust coefficient is determined as follows: 



20 




[Equation 16] 



25 



Where 
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5 



Factor to account for percent change in thrust coefficient per 100°F. This factor is 
0 by default, but can be changed by the user. 

Equation #9 is rewritten to give an expression for thrust, as follows: 



F outi c f out .Pc A^OUti 
[Equation 17] 

10 The variables set forth in the computer model, as shown in FIGs. 6A and 6B are as follows: 
Incremental thrust at output temperature - ThrustOut (i) 
Output temperature - trefout 
Input temperature - trefin 

The subroutine uses equations #10 through #17 to develop arrays of integer numbered 
15 incremental values based on the arrays of values for web burned, burn surface, reference 

pressure, throat area and thrust coefficient, all at the output temperature. The time associated 
with each increment, at the output temperature is derived from the linear web burned and the 
burn rate, at the output temperature and given by the equation: 



f wb outi -wb outi _^ 



[Equation 18] 
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The subroutine solves Equations #10 through #18 for each integer numbered increment 
and also solves Equation #8 for total propellant weight burned. Equation #12 is dependent on 
the pressure, at the output temperature, for the preceding increment, so the process begins with a 
default pressure of 14.7, for the increment / - 1, and writes the pressure derived from Equation 
#15, into each succeeding calculation of Equation #12. The subroutine solves the referenced 
equations, repeatedly, in an iterative operation, varying the burn surface modifier and selecting 
the burn surface modifier which produces a total propellant weight burned, from Equation #8, 
matching the user input propellant weight. The subroutine uses the selected burn surface 
modifier to generate the burnback profile, thrust coefficient versus time, pressure versus time and 
thrust versus time, all for the user selected firing temperature. The burnback profile and thrust 
versus time data are output in graph form, as shown in FIGs. 6A and 6B. Bum surface profiles, 
for various rocket motors, which the method has generated appear as FIGs. 7,9,1 1 and 13. From 
this information, a complete analysis of the ballistics performance of the rocket motor can be 
predicted, for a firing temperature selected by the user, without the burden and expense of test 
firing the rocket motor at the selected temperature. 

The method of the present invention was used to predict the ballistics performance of an 
MK 106 rocket motor using test firing data from a high extreme temperature as the input 
temperature and also to predict the ballistics performance of the same type of motor using test 
firing data from a low extreme temperature. In each case, the output temperature was selected as 
70 degrees. The results are shown in FIG. 8, and represent a close correlation of the two method 
outputs for a 70 degree firing temperature as well as a consistent result for a firing temperature 
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between the two extreme temperatures. The process was repeated for MK 106 rocket motors of 
greater age and produced consistent results, as shown in FIG. 10. FIG. 10 also shows the thrust 
versus time curve of a test firing at 70 degrees, which presents a close match with the prediction 
generated by the method. The same process was performed using MK 1 1 1 rocket motors, with 
results shown in FIG. 12 and on the MK 75 rocket motor, as shown in FIG. 14. 

Having fully described the present invention, it will be appreciated that minor variations 
in the manipulation of the data may be used to produce the same results, without departing from 
the essence of the invention disclosed herein. 
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